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Abstract 

To test kinetic theories, simple and practical setups are proposed. 
It turns out that these setups cannot be treated by Boltzmann's equa- 
tion. An alternative method, called the path-integral approach, is then 
employed and a number of ready-for-verification results are obtained. 

1 Introduction 

One of the fundamental mysteries in statistical physics is that the derivation 
of Boltzmann's equation explicitly resorts to the time reversibility of New- 
tonian mechanics whereas the derived equation seems to offer a plausible 
account for time irreversibility. This logic puzzle has been around us for 
more than a hundred years and a large number of interpretations motivated 
to clear up the issue kept on appearing in the literature. 

Without much attention received, we took another direction and ques- 
tioned the mathematical and physical validity of Boltzmann's equation[lJ. 
In this paper, we concern ourselves with proposing simple and practical se- 
tups, in which Boltzmann's equation, or any other alternatives, can be put 
to computational and/or experimental tests. 

Two types of collisions get involved: particle-to-boundary collisions and 
particle-to-particle collisions. In view of that every particle can be thought 
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of as a particle scattered by boundaries or by other particles, the generality 
of this discussion is rather obvious. 

It turns out that for the proposed setups Boltzmann's equation yields 
either incomputable or unreasonable results. With that in mind, an alterna- 
tive method, called the path-integral approach, is employed and a number 
of ready- for- verification results are derived. 

If interested enough, readers may devise their own ways to prove or 
disprove every conclusion offered by this paper. 

2 Two setups 

To test Boltzmann's theory, as well as its alternatives, we advance the fol- 
lowing two setups, which can be easily realized in real or computer-simulated 
experiments. 

• Setup 1: Referring to Fig. 1, consider a parallel beam whose distri- 
bution is 

f'(v') = n'.g(v' x )5(v' y -0)5(v' z -0), (1) 

where v x > and g{v' x ) is a nonnegative normalized function; and let 
it hit a solid boundary of finite size. It is obvious that the distribution 
given above is quite 'normal'; for instance the relationship 

J f'(v')dv' = J /'(v') dv' x dv' y dv' z = n' (2) 

exists, where n' is the ordinary particle density of the beam. 



j, . . Solid surface 



Paths of 

scattered particles 

Figure 1: Particles of a beam collide with a solid sur- 
face. 

Setup 2: Consider two head-on parallel beams shown in Fig. 2. Beam 
which moves rightwards is described also by the distribution function 
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([T|); and beam 1 is described by 



/IK) = n[ ■ gi (v' lx ) 5(v' ly - 0) 5(v' lz - 0), (3) 

where v[ x < and gi(v' lx ) is again a nonnegative normalized func- 
tion. Though having the same mass, the particles belonging to beam 
and beam 1 (sometimes referred to as type and type 1 particles 
respectively) are considered to be distinguishable in this paper. The 
transverse sections of the two beams are finite, say, shaped like circular 
ones. 



/'— ;;;;; - — n 



Paths of 

scattered particles 

Figure 2: Two beams collide. The collision region is 
highlighted by dots. 

To make the discussion less complicated, we shall ignore any external 
forces and disregard particles that suffer collisions twice or more (letting a 
be small and/or the gas be dilute). With the assumption that all colliding 
parallel beams are produced by stable external sources, the two setups are 
considered to be time-independent. 

3 Application of Boltzmann's equation 

According to the standard theoryOHj, in a spatial region where collisions are 
ignorable, the distribution function, denoted by f(t, r, v) with r representing 
the position (x,y,z), satisfies the collisionless Boltzmann equation 

df df F df 
at or m ov 

in which F stands for the external force acting on the particles. In a region 
where collisions cannot be ignored, the regular Boltzmann equation reads 



5/ +v . df + F df_ 
dt dr m dv 



[ [ 2u(f'f[-ff 1 )a(n c )dn c d-v l , (5) 

J VI Jq c 
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where /'(v') and /{(v^) represent the particles that collide to make /(v) 
increase, /i(vi) represents the particles that collide with /(v) to make /(v) 
decrease, (2u) is the relative speed of the colliding particles , Q c is the solid 
angle formed by the initial relative velocity and final relative velocity, and 
cr(O c ) is the cross section in terms of f2 c . 

In the following two subsections, we shall apply (H|) and © to the setups 
given in Section 2. 



3.1 Application to setup 1 

It is trivial to see how the parallel beam illustrated in Fig. 1 obeys the 
collisionless Boltzmann equation dH). Noting that the beam is stationary 
and there is no external force, namely df'/dt = and F = 0, we obtain 
from (H|) 

, df'{V) df 

v = or — 

or or 



= 0, (6) 

path 



which is indeed satisfied by (JTJ) . For later use, it should be mentioned that 
([6]) characterizes parallel beams. 

We now turn our attention to the particles scattered by the solid surface, 
whose distribution function will be denoted by /(i, r,v). Again, due to the 
existence of df /dt = and F = 0, the collisionless Boltzmann equation (j4]) 
yields 

df 

— = or /|p a th = Constant. (7) 
dr pa th 

The deduction of ([7]) is simple, but difficult things arise. As Fig. 1 has intu- 
itively shown, the scattered particles will not form a parallel beam and their 
density will in general decrease along the paths (which can be easily verified 
in experiments), but according to ([7]) they seem to form a parallel beam 
and their density seems to keep constant along any path of the scattered 
particles. If we wish to adopt d7J, for whatever reason, we shall encounter 
great difficulty in interpreting it, let alone verify it in experiments. 



3.2 Application to setup 2 

In this subsection, /(i, r, v) denotes the distribution function of scattered 
beam particles shown in Fig. 2. 

Firstly, we examine the situation outside the collision region. Again, the 
collisionless Boltzmann equation gives rise to 

v • — = or /[path = Constant. (8) 
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Some of the problems related to ([8]) have been discussed, of which one is that 
the density of the scattered particles will decrease along the paths, while (JSj) 
suggests otherwise. 

Then, we wish to find out what happens inside the collision region, where 
the regular Boltzmann equation (0) is supposed to hold. 

The second term on the right side of ([5]) takes the form 



It is trivial to find that this integral is computable. However, since © is 
associated with particles that suffer collisions two times or more (/ itself 
stands for particles produced by collisions), we shall consider it no more. 
After is omitted, © becomes, by virtue of df/dt = and F = 0, 



Instead of solving (|10p . the following strategy will be taken. Under the 
assumption that / is completely known, say by experimental means, we try 
to find out whether or not (|10p is mathematically computable. Seeming 
peculiar, this strategy serves us quite well. 

The left side can be trivially calculated: at a specified 'position' (r,v), 
the side takes a definite and finite value (/ there is regularly differentiable) . 
The calculation of the right side is nontrivial: to get a definite value at 
the same specified 'position' (r, v) is actually impossible. The basic reason 
behind the incomputability is that v', v^, v, vi and O c are associated with 
each other rather solidly. Note that v has been specified and (v' + v^) points 
along the ±v x axial in this setup; thus vj has no choice but to lie on the 
plane determined by the v x axial and v (under the momentum conservation 
law). Similarly, after v and Vi are specified, there is no degree of freedom 
for Q c . All these tell us that the integration domain of (|10p is just two- 
dimensional and the integral cannot be calculated (See Appendix A, C and 
Ref. 1 for more analytical discussions). 

In fact, there is no great difficulty to configure other types of setups in 
which paradoxes concerning Boltzmann's equation show up. 

4 Application of the path-integral approach 

A sharp and inevitable question arises. What are the true obstacles that 
keep us from mastering the collective behavior of classical particles while 




0) 




(10) 
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each of such particles has been formulated for so long? After a long-time 
effort to find a way out, we come to realize that particles in such system 
involve a dual nature. Before and after they collide, deterministic path- 
equations are obeyed; when they collide with boundaries or other particles, 
indeterministic laws have to be invoked (a is of probability in nature). With 
the motivation to incorporate path-information into an integral of collision, 
a new path-integral approach has been proposed^, [7]. Here, we wish to 
present a simplified version of it. 



Effective cone — AO 

V 




Detector 

v An 

Figure 3: A real or imaginary detector placed at a point 
where the distribution function is of interest. 



The primary concern of the proposed theory is focused on what can be 
measured in experiments. As in Fig. 3, we consider a particle detector placed 
at a position where the distribution function needs to be determined. There 
are key points worth mentioning about the detector. 

1. The detector has a really small opening, whose area is denoted by 
AS . 

2. Without specifying the concrete structure of the detector, it is assumed 
that every scattered particle within a specific velocity domain 

Av = v 2 AvAtt, (11) 

will be registered and any others will not. (The central axis of AO is 
perpendicular to ASq.) 

3. While ASo and Av are considered to be infinitely small, the domain 
of AO, though also small, has been assumed to be finite and fixed. As 
will be seen, the discrimination against AO is taken almost entirely 
from necessity. 
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If the detector works as described above, do we know the distribution 
function at the entry of the detector? The answer is basically a positive 
one. If AN is the number that the detector counts during At, the particle 
density in the phase volume element is 

/('■r...An )a . ( ^ (12) 



where r is the position of the detector opening. The form of (|12p reveals one 
of the most distinctive features of this approach: the distribution function is 
formulated directly. Apart from other merits, this brings a lot of convenience 
to the verification work. 

Examining Fig. 3, we see that only the collisions within the shaded spa- 
tial cone — Afi, called the effective cone hereafter, can directly contribute to 
AiV. Notice that the effective cone is defined with help of path information. 

We shall apply the concepts introduced above to our two setups. 



4.1 Application to setup 1 

Referring to Fig. 1, we wish to determine the distribution function of scat- 
tered particles f(t, r, v, Ail) under the assumption that the detector opening 
is at r and there is a piece of the solid surface inside the effective cone — Af2. 

One of the essential features related to particle-to-boundary collisions is 
that scattered particles will behave like particles emitted by the surface[5]. 
In terms of probabilities, the 'emission' rate p from a surface element dS' , 
whose position is r', can be defined such that 

dN = p(t', r, v', n')dt'dS'dv'dn', (13) 

represents the number of the particles emerging within the solid-angle do- 
main dQ' in the time interval dt' and speed range dv' . Generally speaking, p 
depends on how the particles are 'sent' to the boundary, how many particles 
get involved and what kind of boundary exists there. For purposes of this 
paper, we simply assume that p related to each surface element has been 
known (by experimental means for instance). 

The particles 'emitted' from the surface element dS' at t' will enter the 
detector opening at t and immediately after getting in they occupy 

[r - r'\ 2 AQ, v'dt, (14) 

where the solid-angle domain Af^o is defined as Ailo = AS"o/|r — r'| 2 ; and 
there is a time delay between t' and t 

t = t' + \r-r'\/v'. (15) 



S 



Thus, according to (|12p . the distribution function at r takes the form 



n ,r,V, ) ~JASAvAO \r-r^An v'dt■v^AvAn , U ° j 

where the integration domain AS includes all surface elements enclosed by 
the effective cone —AO (with no blockage assumed). In writing (|16p . we 
have used the fact that AOo is 'much smaller' than AO, namely 

AO < AO, (17) 

which has been ensured by the fact that the detector opening So is infinites- 
imal while AO is finite. With help of (|17p . every particle which moves from 
AS and enters the detector will be considered as one whose velocity is in 
the domain AO. We finally arrive at 

f(t,r,v,An) = -3— / : — . 18 

V s Ail J as r — r r 



It is easy to see that (|18p can be used to calculate the behavior of the 
scattered particles shown in Fig. 1. However, a much more crucial point 
is that (|18|) reveals something that standard approaches have always over- 
looked: discontinuity almost everywhere. When two adjacent effective cones 
pertaining to one spatial point contain different boundaries, the velocity dis- 
tribution there will generally involve discontinuity. In turn, this means that 
any finite-size boundary will generate discontinuity of velocity distribution 
at every path-reachable point (collisions along the paths erase such discon- 
tinuity partly and gradually though). 



4.2 Application to setup 2 

Referring to Fig. 2, we now determine the distribution function of scattered 
beam particles f(t, r, u, AO), where r is located inside or outside the col- 
lision region. 

The collisions taking place within dr' , where r' denotes a spatial point 
inside the effective cone, can be represented by 

[2udr'}[f\t\r'y)dV][f[(t',r'y i )dV 1 }[a(n c )At'}dn c , (19) 

whose derivation can be found in usual textbooks [3j- If the particles pro- 
duced by the collisions expressed by (fT9j) reach r at the time t, there must 
be a time delay expressed by t' = t — |r — r'\/v. 
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With help of the concept of effective cone and the concept of time de- 
lay, path-information has been fully incorporated. Furthermore, the two 
concepts manifest an intriguing relatively new mechanism. Every physical 
event inside the effective cone has its 'right' to contribute to the later-time 
distribution function at r, regardless of its distance to r. With this mecha- 
nism in mind, we shall integrate (|19p over the entire effective cone. (In our 
general theory, the same measure is taken except that each path involves a 
certain loss probability.) 

Since this setup is entirely time-independent, time variables will be made 
implicit in all formulas. 




v 2 AvAQ r 



Figure 4: The relation between the velocity element 
v 2 Av AQq and the velocity element u 2 dudQ c . 

Inserting ([1]) and ([3]) into (fT9|) and then integrating (fT9|) . we finally ob- 
tain, from expression (TPSl) . 



dr'f dn c r d v>[° (20) 

—An JAvA^o Jo J-oo (|r - r'\ 2 AQ, v){v 2 Av Att) 

where AQq = ASq/\t — r'\ 2 and the integration domain (Av, Af^o) implies 
that only the particles in At; and in Af^o can be taken into account. Again, 
since AQq <C Ail holds in most regions inside the effective cone, particles 
that can enter the detector will be considered as ones within AQ. Based on 
the notation c' = (V + v'J/2, c = (v + vi)/2, u' = (V - v'J/2, u = (v - 
vi)/2, d = |c|, d = |c'|, u = |u| and u' = |u'|, the variable transformation 
from (v' x ,v' lx ) to (d,u') and finally to (c,u) is made, and (f20|) becomes 

dv' I dn r r do r d J^ C )\\J\\n'n[ g (c + U ) 9l (c-u) 



-An JAv,An J-oo J\c\ (|r - r'\ 2 Att v)(v 2 AvAtt) 

(21) 

in which the Jacobian between the variable transformation is 

\\J\\ = \\d(v' lx ,v' 2x )/d(d,u')\\=2. (22) 
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By investigating the situation in the velocity space shown in Fig. 4, the 
following relation can be found out: 

/ u 2 dVt c / du ■ ■ ■ « u 2 AvAO • • • • (23) 

JAv,A^o J\ c \ 

Therefore, the distribution function at r is 

where O c is the solid angle formed by u and u' and u is defined as u = 
v(r — r') /|r — r'| — c. 

It is very obvious that (|24|) is truly calculable and the result is a definite 
number. 

Since AO is a finite solid- angle domain, expression (I24p is nothing but 
the distribution function averaged over AO, which is still different from the 
'true and exact distribution function' in the standard theory. However, if 
we let AO tend to zero, AOo <C AO will no longer hold and the concept 
of effective cone will no longer be valid, thus making the whole formalism 
collapse. A thorough inspection of this approach compels us to believe that 
the true and exact distribution function is beyond our reach and any attempt 
to formulate it will ultimately fail. 

5 Conclusions 

The viewpoints and calculations presented in this paper are waiting to be 
proved or disproved by real or computational experiments. If they are right, 
as earnestly expected by us, many fundamental and interesting questions 
arise. For one thing, this paper hints that differential equations should not 
be considered as universally effective apparatuses (as many may assume). 
The following observations are worth mentioning. 

• In a differential formalism, we either respect path-equations, in which 
v = dr/dt holds, or respect partial differential equations, in which t, r 
and v are completely independent. In general, we cannot incorporate 
path-equations into a partial differential approach. 

• In view of that the newly-formulated distribution function is an av- 
eraged one, the proposed theory is 'approximate' in nature and the 
time-irreversibility (information loss) has been built in. 
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• According to the proposed theory, the spatial discontinuity associated 
with a boundary will transform itself into the pervasive velocity dis- 
continuity in the nearby, as well as distant, region. The role that 
such discontinuity plays in a variety of non-equilibrium phenomena 
has been overlooked. 

For the time being, these observations, as well as many related others, 
are just matters for conjecture. Reference papers can be found in the regular 
and e-print literaturejU [7j. 
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Appendix A: Disproof of Boltzmann's equation 

In this appendix, we shall investigate Boltzmann's equation in a brief and 
direct manner (involving no physical assumptions). 

According to the philosophy of standard theory [3], there is an expression 
characterizing how collisions lessen particles around a specific position r and 
a specific velocity v 

«- A Wo : ■ (25) 

At,Ar,Av->0 At ArAv 

where AiV out is the number of the particles that get out of ArAv during 
At because of collisions. Likewise, there is an expression characterizing how 
collisions produce particles around r and v 

At,Ar,Av^0 At ArAv 

where AiVi n is the number of the particles that get in ArAv during At 
because of collisions. Although expressions (f25l) and (|26j) have served as two 
of the basic concepts in our understanding to the statistical world (explicitly 
or implicitly), the recent studies of us fall short of justifying them. 

In fact, a paradox associated with (|25|) can be addressed directly. To 
calculate expression (j25|) . we are supposed to formulate how many particles 
initially within ArAv collide during At and thus leave ArAv. However, 
if the length scale of Ar, denoted as |Ar|, be much smaller than |vAt| (in 
the standard theory there is no constraint on how At, Ar and Av tend to 
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zero), we end up with finding that all the considered particles travel much 
longer than |Ar| during At and almost all the formulated collisions take 
place outside Ar. That is to say, in order for ([25]) to hold its significance, 
we have no choice but to presume |Ar| 3> |vAt|, which is, unfortunately, in 
conflict with Boltzmann's equation itself, whose basic conception is that t, 
r and v are completely independent of each other. 

A similar paradox exists concerning (1260 . To calculate (|26p . we need 
to examine how many particles involve collisions inside Ar during At and 
acquire velocities within Av. However, if |Ar| <C |vAt|, we find that every 
examined particle will, after its collision, escape from Ar in a time much 
shorter than At and only the particles that collide at the very end of At can 
be possibly regarded as ones that get in ArAv during At. 

We now turn attention to whether or not (j26|) is computable (with 
|Ar| 3> |vAt| adopted for the discussion hereafter). Consider that the gas 
of interest consists of many beams and the particles of each beam move 
at the same velocity. Following the standard methodology, our task is to 
examine how many beam i particles in Ar will emerge within the velocity 
volume element Av after colliding with beam j particles. It is easy to see 
that although the scattering velocities Vj and Vj (of beam i and beam j 
particles respectively) have all together six components, they are confined 
to two degrees of freedom since the energy-momentum conservation laws 
serve as four constraints. Namely, in the velocity space, the ending point of 
v,, as well as that of Vj, has to spread over the energy- momentum shell Sij, 
which is of zero thickness and defined by the initial velocities v^ and v'- (of 
beam i and j respectively). Then, a serious problem related to ([26]) is about 
to surface: the number of the beam i particles emerging within Av, labelled 
as (AiV"i n )ij, cannot be proportional to the volume of Av, while (AA^)^- is 
trivially proportional to At and Ar. To this problem, two possible situa- 
tions are relevant. The first is one in which the shell Sy passes through the 
convergence point of Av (to which Av shrinks), and therefore 

hm 4 — I — ir- oc hm — 5 > oo, (27) 

At,Ar,Av^0 AtArAv a^o a 3 

where Av has been assumed to be cube-shaped with side length of a and a 
is the local surface density of beam i particles on Sij . The second is one in 
which S^ does not pass through the convergence point of Av, and therefore 
no beam i particle will emerge within Av as Av — > 0, i.e. 

hm (ff f - 0- (28) 
At,Ar,Av^0 AtArAv 
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Expressions (|27p and ()28|) show that the rate (|26p is not associated with a 
definite value even after |Ar| 3> |vAt| is adopted. (The manipulation of the 
standard theory makes the problem less visible, but not absent. See Sect. 3 
and Appendix C for more.) 

Appendix B: Two types of cross sections 

To understand Boltzmann's equation, it is of considerable necessity to know 
the differences between the cross section in the center-of-mass frame and the 
cross section in the laboratory frame. 

First of all, we examine the collisions between a particle beam (of type 0) 
and a single particle (of type 1) in the center-of-mass frame. Let fi c denote 
the solid angle formed by the velocity of a scattered type particle with 
respect to its initial velocity. Then, the cross section a(Q c ) is defined such 
that 

dN = a(Q c )dQ c [or dS c = a(Q c )dQ c ], (29) 

represents the number of the scattered type particles emerging within the 
solid-angle domain d£l c per unit time and per unit relative flux. Since the 
particles emerging within d£l c actually come from the incident area dS c , the 
cross section has the dimension and units of area. As textbooks have clearly 
elaborated, cr(f2 c ) can be applied normally and easily [2j [3]. 

Now, we examine the cross section in the laboratory frame. In a seem- 
ingly identical way, consider the collisions between a particle beam (of type 
0) with velocity v' and a single particle (of type 1) with velocity v^. Ac- 
cording to textbooks of statistical mechanics, there is a cross section a such 
that [3] 

dN = a(V , v[ — > v, vi)dvc?vi (30) 

represents the number of the scattered type particles emerging between v 
and v + dv, while the type 1 particle emerges between vi and vi + dv\, per 
unit time and per unit relative flux. At first glance, like <r(O c ), a in ([30]) 
can be applied freely. Our studies, however, reveal that a possesses truly 
misleading features. 

As well known, the energy-momentum conservation laws are obeyed in 
a collision. Assuming all colliding particles to have the same mass (simply 
for simplicity) and adopting the notation introduced in Sect. 4, we have 

c = c' = Co and u = u = Uq, (31) 

where Co is the velocity of the center-of-mass and no the speed of particles 
relative to the center-of-mass. For any specific pair (v', v^), Cq stands for 
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three constants and uo for one, and therefore the degrees of freedom for 
(v, vi) are not six but two, as shown in Fig. 5. Namely, in the velocity 
space, the velocities v and vi have no choice but to fall on the spherical 
shell of radius uq, called the energy-momentum shell in this paper. With 
help of this concept, we concluded in Ref. 1 that scattered particles should 
be examined with reference to an area element on the energy-momentum 
shell rather than with reference to a velocity volume element. 




Figure 5: Energy-momentum constraints on scattered 
particles: (a) v' and v' x are the initial velocities, and (b) 
v and V! the final velocities falling on the shell S. 

However, some textbooks state that the cross section a makes sense in 
terms of <5-functions[3]. To see whether it is possible and what is meant by 
it, we rewrite (j30j) as 

dN = advdvi = <r|| J\\dcdu = <r|| J\\dcu 2 dudQ c , (32) 

where ||J|| = \\d(y, v±)/d{c, u)|| = 8 is the Jacobian between the two sys- 
tems. The comparison between (j32[) and (|29p yields 

<t(v', v[ -> v, vi ) = 1 a(0, c )5 3 (c - c )5(u - u ), (33) 

u IM II 

where <5 3 is the symbol of three-dimensional <5-function. Obviously, the 
energy-momentum conservation laws have been included in (|33|) . 

Can the cross section a so defined offer extra advantages? Or, on second 
thoughts, can it still cause ill treatments in certain contexts? To answer 
these questions, we look at Fig. 6, in which two small domains for scattered 
particles are specified as Av and Avi (initial beams with v' and v[ have not 
been depicted). Under the assumptions that Av and Avi are symmetric 
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with respect to the center-of-mass and Av encloses a small part of the 
energy-momentum shell, whose area is denoted by AS, we obtain, by virtue 
of {32]) and {H, 

/ ■■■&dvdv 1 = [ ■■■a(Q c )dn c = [ ■■■a{£l c )dS/u 2 , (34) 

JAv,Av! JAn c JAS 

where Afi c is associated with Av, via An c = AS/u 2 . Equation §M§ illus- 
trates that a is essentially defined on the energy-momentum shell S and by 
the cross section <t(S7 c ). 




Figure 6: The velocity domains Av and Avi have to 
be symmetric with respect to the center of the energy- 
momentum shell. 

If there exists no explicit constraint on v and vi (other than the energy- 
momentum conservation laws), expression (|34p becomes 

/ • • • odvdvi = ■ ■ ■ a(Q c )dft c , (35) 

Jv,vi J (An) 

where the subindex (47r) of the last integral means that the integral is over 
the entire domain of Vt c . 

If Av encloses a small part of the energy-momentum shell AS" = u 2 AVL c 
while vi is subject to no explicit constraint, then the expression 

• • • advdvi (36) 

Av,vi 

seems to differ from ()34j) in an obvious way. However, the difference is 
a misleading one. Due to the existence of energy-momentum conservation 
laws, the confinements of v and vi are intrinsically connected, and therefore 
(|36j) yields exactly the same result as (j34"l) does. 
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In Appendix C, it will be unveiled that, associated with mistreating (|36p. 
the standard formalism involves a mathematical error. 



Appendix C: An error in the textbook treatment 

According to textbooks, the regular Boltzmann equation is based on [3] 

?L dl + F,d£ = Hm (AAQ in - (AAQ out 
dt V dr m dv AtArTv^o AtArAv 

where (AjV); n and (A-/V) out stand for the particles that enter and leave 
ArAv during At due to collisions. In this appendix, we ignore the paradoxes 
unveiled in Appendix A and concern ourselves solely with how textbooks 
treat (|37p in the mathematical sense. 

The textbook methodology of formulating (AA r ) out is that, for a defi- 
nite Ar, the particles leaving Av during At due to collisions are identified 
as (AiV) out . On this understanding, /i(vi)dviAr represents the particles 
knocking some particles out of the beam f(v)dv, and we get 

(AA0out=/ 2uAr[f(v)dv][f 1 (v 1 )dv 1 }[&At}dv'dv' 1 , (38) 

where a = <r(v,vi — > v',v^) is the cross section in the laboratory frame. 
With help of (f35l) . we get 



sJT n itf A"' = / 2na(O c )/(v)/ 1 (v 1 )dv 1 dfi C) (39) 

AtArAv->0 At ArAv ,/vi ,0.(470 



/vi,Oj(4ti-) 



where vi is unlimited and the domain of f2 c is (47r). So far, every thing is 
exactly the same as that in the textbook treatment. 

A similar methodology is employed to formulate (AiV)j n . f'(v')dv' and 
/{(v'^dv^ are identified as two colliding beams, and we have 

(AN) in = f 2 U Ar[/ , (v / )dv , ][/((v / 1 )dv' 1 ][ ( rAt] ( ivdv 1 , (40) 

J v'jVj , Av,vi 

where a = <r(v',v^ — > v,vi) and Av is used to remind us that only the 
particles emerging within Av will be taken into account. Although the 
symmetry between (|38l) and (|40l) appears obvious and has played a vital 
role in deriving the Boltzmann equation, we are now convinced that it is 
an illusive one. Due to the existence of the energy-momentum conservation 
laws, when v in (140 p is limited to an infinitesimal volume Av, vi must be 
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limited to the symmetric Avx, while vi in (|39p is unlimited. That is to say, 
(|4"D"]) should be replaced by 



(AN) in = f [2«Ar][/'(v / )dv / ][/{(v / 1 )«ivi][ffA*]dvdvi. (41) 

Jv',v^,Av,Avi 

So, the aforementioned symmetry disappears and we arrive at 

. A n ™ „ + I 2ua(n c )f'(V)f[(V 1 )dv 1( m c . (42) 

AtArAv-^o At ArAv JvinJin) 



As Appendix A has shown, a serious evaluation of the left side of (|42p 
actually leads us to nowhere. 
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